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Abstract: In this work we present a systematic construction of the potentially ghost- 
free non- linear massive gravity actions. The most general action can be regarded as a 
2-parameter deformation of a minimal massive action. Further extensions vanish in 4 
dimensions. The general mass term is constructed in terms of a "deformed" determinant 
from which this property can clearly be seen. In addition, our formulation identifies non- 
dynamical terms that appear in previous constructions and which do not contribute to the 
equations of motion. We elaborate on the formal structure of these theories as well as some 
of their implications. 
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1. Introduction and Summary of Results 

There are several motivations for studying massive gravity. From the observational point 
of view, the conclusion that 95% of the universe is made of dark energy and dark matter 
is closely related to the assumption that general relativity is equally valid at all length 
scales. It is possible that modifications of general relativity over large distances change 
this picture and massive gravity is one such modification. It is also possible that a graviton 
mass ameliorates the cosmological constant problem. From the theoretical perspective, the 
shear difficulty of constructing a consistent theory of massive gravity makes the problem 
more interesting. 

In 1939, Fierz and Pauli |l], ||] constructed the ghost-free theory of non-interacting 
massive gravitons in a flat background. Since then, this theory has stubbornly resisted a 
non-linear generalization. The work of Boulware and Deser || showed that a generic exten- 
sion of the Fierz-Pauli (FP) theory to curved backgrounds will contain ghost instabilities. 
The ghost analysis in massive gravity was greatly simplified by the work of Q, motivated 
by the Goldstone-vector boson equivalence theorem in quantum field theory. The authors 
in [|[] argued that, since a graviton mass can be associated with a spontaneous breaking 
of general covariance, the associated Goldstone modes in a non-unitary gauge will contain 
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the ghost information. Based on this observation, a procedure was outlined in [Q, |5| where 
the ghost could be avoided order-by-order by tuning the coefficients in an expansion of 
the mass term in powers of the metric perturbation hnv and the longitudinal Goldstone 
mode 7r. An added advantage is that within this effective field theory-like framework, if a 
ghost turns out to be heavier than the UV cutoff of the effective theory it can be safely 
ignored. Thus even when the ghost is not completely avoidable, it could still be manage- 
able. The preliminary conclusion in was that at the quartic order in Goldstone modes, 
the expansion will unavoidably contain ghosts. 

This problem was recently revisited by de Rham and Gabadadze Q who showed that 
massive gravity ghosts can indeed be avoided, or at least their appearance relegated to 
higher energy scales where the effective field theory may no longer be applicable. One can 
then successfully construct a theory that is ghost-free in the so-called decoupling limit, 
which amounts to retaining terms only first order in h^ u , but all orders in ir || 0. In 
particular, Q constructed a non-linear massive action starting from the decoupling limit 
data. Using this theory it becomes possible to study the ghost away from the decoupling 
limit. It was concluded that one of the resummed theories is ghost-free at least up to 
quartic order in the nonlinearities. 

In this work we present a systematic construction of massive gravity actions that 
simplifies the analysis and brings out several important features. The main results are 
summarized below. The implications of these results are discussed in more detail in the 
concluding section. 

• We construct the most general, potentially ghost-free, mass term as a "deformed" 
determinant. The general massive gravity action can be regarded as a 2-parameter 
deformation of a minimal massive action which is, in turn, the minimal non-linear 
extension of Fierz-Pauli massive gravity. The origin of the 2 parameters can be easily 
seen in this formulation. No further generalizations exist in 4 dimensions. The full 
nonlinear equations of motion are also written down. 

• It is shown that the earlier constructions of non-linear massive actions lead to non- 
dynamical terms in the action that do not contribute to the equations of motion. 
These terms arise due to an assumption made to simplify the construction. In our 
formulation it becomes straightforward to isolate and remove these terms, implying 
that the resummed theory found in Q has an equivalent description in terms of a 
simpler action. This also shows that while the massive action was obtained starting 
with a covariant version of the Fierz-Pauli mass, it also turns out to describe theories 
which reduce to the non-covariant Fierz-Pauli mass at lowest order and, hence, is 
more general than it seems. 

• We emphasize the physical relevance of generalizing the auxiliary metric r]^ u , used 
in the construction, to a general f^ u . In this context, we discuss solutions of the 
minimal theory that have the potential to avoid the observational consequences of 
the vDVZ discontinuity ||, ^] as well as screen a large cosmological constant. 
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The actions presented here are known to be ghost-free at the interactive level to first 
order in the metric perturbation hn V and to all orders in the Goldstone longitudinal mode, 
7T. This was extended to quartic order by an ADM analysis in 0. However, a proof that 
these actions are ghost free at the complete non-linear level is still lacking (see, for example, 
|1C] for an argument that the ghost is in fact present at the quartic order). But the stronger 
statement is that if any such non-linear ghost free action exists, it will have to belong to 
this family of actions (see note added at the end). 

The organization of this paper is as follows. In section 2 we review the basic features 
of massive gravity, paying particular attention to the features that will be most relevant 
when constructing a nonlinear theory. We also briefly summarize the recent developments 
in constructing ghost-free nonlinear extensions of the FP action. In section 3 we construct 
the minimal non-linear extension of the Fierz-Pauli mass term and discuss its screening 
solutions as well as the origin of non-dynamical terms. In section 4 the generalization of 
this model is considered. We describe the formal structure of the theory and argue that 
there exists at most a 2-parameter family of actions with the desired properties. We write 
down the full non-linear equations of motion for the generalized theory. Some implications 
of our results are discussed in section 5. 



2. General Structure of non-Linear Massive Gravity 

In this section we review features of massive gravity that are relevant to the construction of 
non-linear mass terms, and make some observations that clarify the underlying structures. 
We then briefly discuss recent developments in constructing ghost-free non-linear extensions 
of Fierz-Pauli massive gravity. This provides the background and conventions for the 
systematic approach of the following sections. 

2.1 The Fierz-Pauli mass in linear theory 

Fierz and Pauli |l|, constructed a mass term for linearized metric fluctuations around 
flat spacetime, = g^ u - r/^, 

9 

— (v^" - a K h ») ' with a = 1 - ( 2 - x ) 

For any other value of a the theory has a ghost. Invariance under infinitesimal general 
coordinate transformations (GCT) can be restored by introducing Stiickelberg fields and 
replacing h^ u by h^ u + d^ir^ + d^^. In the Stiickelberg formulation, the longitudinal mode 
7T defined through 7r M = ir^ + d^ir, appears in the action through a 4-derivative term and 
carries the ghost degree of freedom. The ghost decouples only for a = 1. Unfortunately, 
unlike the case of massless linearized general relativity, there is no obvious unique way to 
extend the massive linearized theory to the non-linear level. But in recent work, ghost 
decoupling has turned out to be a powerful criterion. Let us first discuss the general 
features of a non-linear mass term before discussing the recent developments towards its 
construction. 
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2.2 The auxiliary metric and the inverse metric 

On general grounds, a non-linear mass term in gravity (a) cannot contain derivatives of g^ u 
and (b) should be invariant under general coordinate transformations (GCT). Such a term 
cannot be constructed using the metric alone since the only two relevant quantities, det g 
and Tr<7 = 4, are not adequate for this purpose ||. To overcome this, one of the following 
options may be used: 

1. Introduce an extra rank two tensor f^ u as in 0,0]. 

2. Introduce an auxiliary extra dimension |11, 12]. 

3. Generate a mass through a gravitational Higgs mechanism 1 [|l3[ 14, [D| [l(], [D], [Hf . 

Both options (2) and (3) can be formulated in terms of (1) as will be discussed below for 
Higgs gravity and in []l9| for the extra auxiliary dimension model. Hence we will concentrate 
on this first case. 

Combining the dynamical metric g^ v with a covariant rank two tensor f^u(x) (hence- 
forth called the auxiliary metric), one can construct arbitrary scalar functions F{g^ 1 f) 
and F'(f~ 1 g). The GCT invariant mass term has to be written in terms such functions. 
A special case is when the mass term involves only F{g _1 f), i.e., it depends only on the 
inverse metric g^ v . A simple inspection shows that the massive gravity actions considered 
in [||], and in the subsequent works || ^, [?]], all belong to this class 2 , even though, in 
practice, the setups may have been described in terms of g^ v . In the last paragraph of this 
section we argue that the combination g' 1 f is adequate to describe any ghost free mass 
term and so will be explicitly assumed here. Thus we consider massive gravity actions of 
the form, 

S = -M 2 J d A x^R{g) + M 2 p m 2 j d 4 Xy/^g F(g~ 1 f) . (2.2) 

Invariance under general coordinate transformations implies that F must be expandable 
as a sum of products of the type Tr[(g _1 /) n ]. 

Invariance under GCT is insured since transforms as a rank 2 tensor. To make 
this explicit, one can parametrize the auxiliary metric as a coordinate transformation of 
some fixed metric j|], 

d(b a - d6 b 

M*) = £ £ ■ ( 2 -3) 

The (j) a (x) are to be regarded as general coordinate transformations of the x^ and generalize 
Stiickelberg fields to non-linear level. f a b is a fixed auxiliary metric that should be chosen 



1 Note that in this context, the term Higgs mechanism does not imply the existence of a mass generation 
mechanism on par with that in vector gauge theories, and the corresponding ramifications for unitarity, etc. 
Here, it mainly refers to the association of the graviton mass with scalars of the form </> M = x M + 8cj>. We 
retain the terminology with this caveat. 

2 This also holds for the "auxiliary extra dimension" model of jll], [Til which is special in that it 
admits a description in terms of both g -1 / and WM- 
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consistent with the physics of the problem. Both these are scalar functions 3 of (/ is a 
tensor with respect to reparametrizations of 4> a ). The physical (unitary) gauge corresponds 
to choosing cp a (x) = x a . In ||] f is identified with a background dynamical metric g in the 
physical gauge. If / is taken to be a flat metric, the cp a can be chosen such that f a b = r]^. 

2.3 Non-linear generalizations of the Fierz-Pauli mass 

Let us now describe the ideas that have recently been used to construct non-linear graviton 
mass terms 4 . Given the auxiliary metric f^, a covariant non- linear version of the Fierz- 
Pauli mass (^4]) can be constructed in terms of a matrix H„ given by, 

9^f^ = SH-HS}. (2.4) 

This is obviously a (1, 1) tensor under GCT. Note that, rewritten as H^ v = g^xH^ = 
9fj.v — ffiu, one recovers the matrix introduced in [Q] in a more elaborate setup. Any 
mass term can be expressed entirely in terms of H„, and we will only retain this in the 
description. Then the Fierz-Pauli mass can be generalized to a covariant non-linear form 
(see § and references therein), 

2 

' (nam -mil*), (2.5) 



4 

which reduces to the original FP mass (|2.1| ) at the quadratic order in h^ v and in the unitary 
gauge for <p a . We refer to this as the covariant FP term. 



But this is not enough. Though (2.5) is ghost free at the quadratic order around flat 
background, the ghost mode generally reappears at the quadratic order around curved 
backgrounds p]. Equivalently stated, the ghost reappears at higher orders around flat 
background. To cure this problem, the idea, as proposed in EJ and pursued in ||, ||], is to 
add terms to (^1^) that are higher order in Hy to cancel the ghost at higher orders. 

To simplify the ghost analysis, the theory is often written as an expansion around flat 
background, which for g 1 *" gives, 

g ^ = + . (2.6) 

In terms of h^ v = g^ v — rj^, one has f)'"' = Yl^=i [( — r]~ 1 h) n " r l~ 1 }^ U ■ Also, since the <p a {x) 
are general coordinate transformations of x^, they can be expanded around the "unit 
transformation" as pj, 

(p{x) = a" + 7r"(x) . (2.7) 

The 71"^ transform as "coordinate differences" under GCT. Then Hit given by Q2.4j ) take 
the form (on choosing f^ v = to facilitate comparison with FP theory), 

Hi = -3*v„ - d v ^ - d^ a d U TT a - w {8p + d u ir p + d p 7T U + ^vr a cu a ) (2.8) 



3 To see this, start with a coordinate system {<^ M } and consider two reparametrizations {a;''} and {x 1 ^} 
defined by M (a;) and 4>'^(x'). Now, a coordinate transformation from x to x' is subject to M (x) = 
clarifying that ^(x) are scalars. 

4 While the structures introduced in this subsection are useful for motivating the form of the non- 
linear actions of later sections, and for seeing the relation with other recent work, they are not needed for 
formulating the non-linear theory, which is good news in view of footnote H 
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Here, all indices have been raised or lowered with respect to their natural positions with 
the background metric rj^. The expression is exact in tt^ and in \f" v which, so far, need 
not be infinitesimal. If the graviton mass is assumed to arise as a result of a spontaneous 
breaking of GCT, through some kind of gravitational Higgs mechanism, then the TT a would 
be the unphysical (pure gauge) Goldstone modes associated with the symmetry breaking 
0,0- 



The field most relevant to the ghost problem is the "longitudinal" mode tt, given by, 

= ^ + . 

Let's use the notation OH for the matrix with elements Hy to avoid mixup with H = H^. 
Also, let II denote the matrix with elements Yiy = d^d v ix. Then retaining only tt, and 
setting tt^ = 0, one has, 

g- 1 f = (l-h) = (t + t }V )(t + U) 2 . (2.9) 

This expression is exact 5 to all orders in t) and tt. 

Recent developments: As discussed in [|, the Boulware-Deser ghost instability of 
massive gravity is associated with the appearance of tt dependent terms in the action with 



more than 2 time derivatives. Thus, in the massive gravity action (2.2), an expansion of 
F(l — bi) in powers of OH should be constructed in such a way as to avoid this problem. Such 
a program was carried out successfully by de Rham and Gabadadze in ||. They obtained 
a 2 parameter family of actions to quintic order in OH that avoids the ghost problem at 
least up to first order in h au and all orders in tt. Subsequently Jjj proposed a method for 
resumming such expansions in powers of OH into closed-form non-linear actions. 

In the next two sections, we present an alternative formulation of the massive action 
based on a "deformed determinant" that clearly brings out the systematics, showing in 
particular that the framework is more general than might appear at first sight. 

Adequacy of the description in terms of g^ 1 f: Consider a mass term that also 
depends on f g = (1 — H) _1 . It is obvious that any scalar function F'(f~ 1 g) can also be 
expanded in powers of OH. The quintic order polynomials of OH analyzed in || are general 
enough that they do not exclude F' in the mass term. However, as can be seem from 
this work as well as from [f?j, all ghost free combinations are obtainable from actions that 
contain F(g^ 1 f) alone, owing to the simple relation between n and OH in ([2.9|), Hence 
this structure is adequate for describing ghost free massive gravity actions. Though, up 
to a given order, the H-expansions may also be obtainable from an F'(f~ 1 g), that would 
correspond to a complicated rewriting of F(g _1 f) to that order. 



J A word caution on covariance issues: The above parametrization of fluctuations refers to a perturbative 
expansion around flat background. At the non-linear level the -k^ are coordinate differences and hence tt^ 
and 7r are not proper vectors and scalars. Thus, strictly speaking, statements that hold for vector and 
scalar modes at the linearized level, may not fully apply here. Also, the resulting expressions cannot be 
covariantized by simply using covariant derivatives. If one insists on using this type of parametrization 
for the non-linear theory, while maintaining manifest general covariance, then an option is to use Riemann 
normal coordinates and express 7r M in terms of a geodesic tangent vector tt M and derivatives of the curvature 
tensor as in[pl[. 
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Relation to Higgs gravity: In Higgs gravity setups [15, 16, 17, 18], one considers 
gravity coupled to 4 scalar fields, say (j) a with a Minkowskian internal metric. The scalars 
are supposed to acquire a vacuum expectation value < <f) a >= x a through some dynamics 
and break the GCT symmetry of the action, rendering the gravitons massive 6 . Note that 
all terms involving derivatives of the <p a can be written in terms of g~ l f ', or using the trace 
structure, in terms of d /1 (j) a g^ v d v (j) h 'fj, c . For example, the basic kinetic term is Tr(g -1 /). 
This accounts for the practical similarity between Higgs gravity and massive gravity, as, 
for example, in [p3|| . 

3. The Minimal non-Linear Massive Gravity Action 

In this section we construct a minimal non-linear extension of the Fierz-Pauli mass term 
and discuss some of its main features. This is useful for understanding the more general 
non-linear extension to be given in the next section. 

3.1 The mass term as a minimally deformed determinant 

The construction of the non-linear graviton mass initiated in |4|, |5| and successfully em- 
ployed in H [7j , is based on two simple criteria. 

1. The absence of ghosts at low orders: The function -F(<? _1 /) in the massive grav- 
ity action ( |2.2| ) is constrained by the requirement that its expansion around a flat 
background leads to a ghost-free theory. The expansion involves \f" v and ir^ = d^ir 



(ignoring the "vector" modes) through (2.9). As explained in M, the ghost is associ- 



ated with 7r terms with higher than two derivatives which, therefore, should not arise 
in the expansion of F. In particular, the ^-independent part, F(g~ 1 f )\h=o, should 
contain the potentially dangerous ir terms only through total derivatives. The next 
step is to check that the term linear in f) is also ghost free (this is, of course, necessary 
but not sufficient for avoiding the ghost to higher orders in f)). 



2. The recovery of the covariant Fierz-Pauli action ( |2.5| ) at the lowest order. 
Note that the second criterion demands more than the minimum requirement of repro- 



ducing the original FP mass in flat background (2.1). Hence it could, potentially, exclude 
non- linear actions that reduce to ( |2.1| ) without producing ( |2,5[ ). We will show that these 
two classes of actions only differ by non-dynamical terms that can be easily isolated. 

In & ^, 0, to implement the first criterion, total derivative functions of H„ = d^d v Tr 
are explicitly constructed in terms of traces of powers of n ( |2.9| ). In contrast, the approach 
here is based on the observation that all such total derivative terms naturally appear in 
the expansion of det(l + n). This simplifies both the procedure and the final results by 
identifying the redundancies that result from the second criterion. 



6 A dynamical realization of Higgs gravity is provided by the "brane induced gravity" construction, as 
discussed in p3] . While the obvious breaking of translation invariance in directions transverse to the brane 
generates massless Goldstone modes, a more subtle spontaneous breaking of GCT along the brane occurs 
as a result of integrating out bulk modes. This automatically generates the cj> a and the graviton mass in a 
gauge invariant set up, though a ghost free construction is still lacking. 
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To satisfy the first criterion, consider F(g 1 /) at lowest order, i.e., for f) = 0, when, 

GrV)Lo = (i + n) 2 . (3.1) 

As mentioned above, at this order tt must appear in the Lagrangian only in total derivative 
terms. All such total derivative terms appear automatically (with given coefficients) in an 
expansion of the determinant, 

4 _i 

det(l + II) = J2 n , (4 _ n) , ^-^a„ +1 ...a 4 e^ x ^-^d v ^ ■ ■ ■ d^d Un rr . (3.2) 
n=0 

The total derivative nature is manifest due to the antisymmetry property of the epsilon 
tensor. We can choose F\, „ to have this form, but the first two terms of the summation, 

lh=0 ' ' 

1 + d^d^ir, (for n = and 1, respectively) cannot arise in an action that, by the second 
criterion, reduces to the covariant FP form (|2.5| ) at the lowest order (since this is at least 
quadratic in tt). Thus choosing F\^ = q consistent with the second criterion amounts to 
deforming det(l + II) by removing the first two terms from the above expansion, 

F(g- l f)\h=o ~ det(l + II) - Tr(l + IT) + 3 . (3.3) 

This is the minimal deformation of the determinant structure. One may also consider 
non-minimal deformations by changing the coefficients of the n > 2 terms in ( |3.2| ) which, 
obviously, still preserve their total derivative structure. Here we concentrate on the minimal 
deformation, leaving the general case to the next section. 

Now, the field tt enters the action only through g -1 /, in the combination fl2.9|) . In 
particular, det(l + II) = \/det{g~^f)\ h=Q and Tr(l + II) = Tr y/ g~ 1 f\ h _ Q . Thus, going 
beyond h = requires making the substitution 1 + II — > \J g^ 1 f leading to a massive 
gravity action (|2.2j), with the non- linear mass term 



m 2 j d A xyf=gF = 2m 2 j d A x^g yj&et{g- l f) - Tr vV 1 / + 3 • ( 3 - 4 ) 

The construction involves the square root matrix defined such that 7 \[E^fE = E. Again, 
the mass term is given by a minimal deformation of det(y / <7 -1 /) = det(l + \J g _1 f — 1) 
that involves removing the first two terms from its expansion. 

To confirm that our massive action is the non-linear completion of one of the ghost 
free quintic order polynomial actions of ||, we expand it in powers of OH. To lowest order 
in IH, one recovers the covariant form of the FP mass term (|2.5|) . To quintic order in OH 
one obtains the action presented in || with the particular coefficients 8 C3 = | , c% = — ^ 
(fj there multiplies terms that vanish in 4-dimensions as will be discussed in section 4). 
Generalizations of this action are considered in the next section. 

The above action can be trimmed further: 



7 For the metric, we stick to the conventional notation y/—g = y/~ det g. In all other cases, y/E represents 
the square-root matrix and not the determinant. 

8 To correctly match the coefficients at quintic order we note that there is a typo in [H|. The expression 
for /3 should be corrected to /3 = |c3 — 3ds + 20./V. 
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1. The term is needed to cancel out a cosmological constant of order M 2 m 2 . If 
this cancellation is not required, the factor 3 can be replaced by a free parameter. 

2. The first term simply reduces to y '— det f and does not contribute to the metric 
equation of motion. As shown in the next subsection, it also does not contribute 
to the (ft a equations since these are already contained in the metric equation. Thus, 
being non-dynamical, it can be dropped from the action at the classical level (as long 
as one does not intend to promote f^ u to a dynamical variable). The origin of this 
term will be explained below. 

This leaves us with an action that turns out to be the minimal non-linear extension of the 
massive Fierz-Pauli action in the presence of a cosmological constant A = A' + 6m 2 , 



<? • — — /If 2 



'I J d 4 xy^g R + 2m 2 Tr yV 1 / + A' . (3.5) 



Note that although the action is constructed as a minimal deformation of a determinant, 
the determinant itself becomes non-dynamical and only the "deformation" stays on. This 
action is also "minimal" in the sense of having the simplest structure as compared to the 
generalizations constructed in the next section. 

If we now expand the mass term Q3.5| ) in powers of IH we find, to quadratic order, 

2m 2 (1 - \ Tr H - \ Tr(H 2 ) + 0(H 3 )) , (3.6) 

which is not of the covariant FP form (^B|). However, expanding this in the metric per- 
turbation around flat spacetime does reproduce the basic FP form ([O]). Hence, this 
is a valid mass term. What has happened is as follows: As stated above, the first term 
\J — det / in (3.4) does not contribute to the equations of motion. However, this term is 



needed to reproduce the expansion that starts with the covariant FP form Q2.5| ) at the 
quadratic level, as was required by the construction. But this requirement can be relaxed 
since it is totally acceptable that a non-linear massive action directly reduces to the basic 
Fierz-Pauli form (|2.1| ) at the quadratic level in h^ v without ever reproducing (|2.5|). Thus, 



one concludes that while taking the covariant FP term (2J3) as the starting point simpli- 
fies the construction, it does not make it less general. Rather, it only leads to redundant 
non-dynamical terms in the action that are present for "cosmetic" reasons. The formula- 
tion presented here makes this manifest and conveniently isolates the redundant terms. To 
drive home the point, note that in a formulation that does not take care of the redundancy, 



the same minimal theory looks much more complicated than (3.5), i.e., it is given by i 
where \J det(g -1 /) is expanded in traces of powers of the square-root matrix -y/g -1 /. 
Returning to the ghost issue, when expanded in powers of f) around flat background 



the lowest order term in (3.4) or (3.5) is ghost free by construction. The next term, linear 



in f) and to all orders in 7r, is also healthy having the form (see footnote ^| below), 

l -m 2 ^ v (d^-n^U)^. (3.7) 
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In particular no higher powers and derivatives of tt arise at the linear order in f). This 
corresponds to one of the models constructed in Q to quintic order in DH = 1 — <? _1 / '. It 
was pointed out there that, in the absence of non-linear tt interactions, the theory will 



not exhibit the Vainshtein mechanism [24, 25] which is associated with the strong coupling 



behavior of tt in the limit m — > 0. More precisely, non-linear tt terms at higher orders in F) 
may still generate a strong coupling scale but this will be suppressed by the Planck scale, 
leading to a much smaller Vainshtein radius. The issue of ghosts at higher orders in f) is 
still an open question and will not be addressed here. 

3.2 Equations of motion 

In this subsection we write down the g^ u and 4> a equations of motion and, in particular, 



show that the "cosmetic" yj— det / term in (3.4) does not contribute to the 4> a equations. 
The metric equation: This can be worked out using ( |3.5| ) and reads 9 , 

Rftv - ^9fiuR - \^9^v + \m 2 \g\fg- 1 ! + {yfg~ r f) T g 

L L L L J fxu 



+ m V 3 " Tr Vg^f) = G N T^ V . (3.8) 



At least for a diagonal metric ansatz, as in [26], the square-root matrix simplifies and one 
gets closed-form non-linear differential equations that could in principle be solved. 

The left hand side is not automatically divergence free, rather its divergence yields, 



V^vTVlU/VV^vV 1 /); - 2^Tr( v V T 7) = 0. (3.9) 



Here the covariant derivatives are to be written out regarding \J g~ x j as a (1, 1) tensor, 
similar to g^ 1 f ■ This follows from the expansion of the radical in powers of DH, or from the 
fact that, using E = (v 7 ^) 2 , a transformation E — > AEA^ 1 implies y/E — > A^/EA~ 1 . The 
trace of ( |3.8| ) gives, 

R + 2A + 3m 2 (Tr \fg^ r j - 4^ = G N T£ . (3.10) 

A simple class of solutions will be considered below. It is instructive to contrast these 
equations with their linearized forms around flat spacetime (the FP case). The linear 
version of (|3.9j ) implies the vanishing of the linearized curvature scalar Rl . A non-zero Rl 
cannot be recovered in the limit m 2 — > which is at the origin of the vDVZ discontinuity 
of linearized massive gravity. 

The <p a equations: These are equivalent to the divergence of the metric equation 
Q3.9|) . This is a consequence of the general covariance of the action S m i n [g, cfr]. Indeed, under 
a general coordinate transformation, Sx^ = 5g^ u = — 2V^^ and 8cj) a = —^d^ 11 . The 
invariance of the action then implies J d 4 Xy/^g[(-^-^ r7 )V f1 ^ 1 ' + -^=-^d u (j) a ^ u } = 0, or 
on integrating by parts, 

V^(^L^) = ^L#5^. (3.11) 



9 Note that E = (Ve) 2 implies 5E = (8Ve)V~E + \/E(8VE), so that Tt(8VE) = (l/2)Tr{(y r E)- 1 SE]. 
This last equation holds only under the trace. Then for E = g~ l f , a variation of the metric gives 
T^(5\Jg~ x f) — (1/4) Tr [fff\/ 9~ x f + (\/5 _1 /) T ffj <5<? _1 j , where use is made of the cyclic property of the 
trace and /(-y/g -1 /) - = g(g~ /)(\/ff -1 /) - = ff\/<? -1 /- The result is then explicitly symmetrized. 
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The left hand side is the divergence of ( |3.8[) . Since (p a (x) are non-singular coordinate 
transforms of x 1 *, the matrix d v <$> a is invertible. Hence the (fi a equations 5S/5(p a = are 
equivalent to the divergence equations ([?!]). As a consequence, the cj) a equations are not 
affected by \J — det / in (3.4), since it does not contribute to the metric equation. It is 



instructive to check this explicitly. The (p a equations of ( |3.4|) , including v 7 — det /, are, 

V^(D^d u rfac) + D^d^j^-d^ = , (3.12) 

where, 

= Vdet^- 1 /)!/- 1 )^ - \ IV^f + g-WtrHfgY (rV- (3.13) 

l L J p 

Multiplying by the invertible matrix <9\</> c an d manipulating, this reduces to the equivalent 
expression ( |3.9D , as expected. In particular, the contribution from the first term of D^ v 
completely drops out, once again showing that y/— det / is non-dynamical. 

3.3 The screening solutions 

For the massive gravity actions considered here, the proof of absence of ghosts at low 
orders is valid only for a flat auxiliary metric, /«„ = rju U . This theory exhibits the vDVZ 
discontinuity ||, ||] at the linear order which implies radically different bending of light near 
massive objects as compared to general relativity. This difference needs to be avoided at the 
non-linear level, at least in the vicinity of massive objects, should a theory of light gravitons 
describe the real world. In general, massive gravity can achieve this via the Vainshtein 
mechanism 10 [24, 25). However, the analysis of [||, performed in the so-called decoupling 
limit approximation, indicates that this mechanism will not work precisely for the minimal 



model (|3.5|) , as briefly discussed following equation (|3.7|) . Besides this, the divergence 
conditions ( |3.9|) are constraining enough that, when working in a physical gauge, they could 
rule out spacetimes of physical interest, including homogeneous, isotropic FRW cosmologies 



1 28 ] . Some of these issues can already be seen in |29], 30]. Thus, the classical solutions in 
massive gravity models with flat may drastically differ from general relativity. 

Of course, the action ( |3.5D can also be used for non-flat /, although in this case, the 
ghost problem has to be investigated afresh. Let us split the metric into background and 
fluctuation parts, = g^ + and choose the physical gauge, ty* = x M . Now take, 

Uu = c 2 g^. (3.14) 

The caveat here is that, classically, there is no real distinction between the "background" 
and the "fluctuation" (unless T^ v has a natural split in this way). In principle, this problem 
can be avoided by giving appropriate dynamics to f^ v . We will not attempt that here, 



10 But even in that case, the vDVZ discontinuity can be tested far away from massive sources. For example, 
by comparing mass estimates for galaxies and galaxy clusters as derived from gravitational lensing, with 
dynamical mass estimates, such models can be constrained observationally |E7j. 
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except to say that in that case the \J— det / term can no longer be disregarded. If we 
simply add a term M 2 m 2 Xf\/— det / to the action, then the minimization condition, 

(VfT 1 /)/- 1 + rHVg^ff - 2A / Vdet 5 -i/r 1 = (3.15) 
has the solution c = (1/Aj)s. More generally, the above ansatz naturally arises as a 



solution in bi-gravity theories |yj 32|. Here, we will not attempt giving any specific 
dynamics, but will bear this picture in mind in what follows. The intention is to discuss 
the properties of solutions that resemble ( 3.14| ) in the appropriate extension of the theory. 
With the above ansatz, the divergence condition ( [3.9| ) gives d^c = and the metric 



equation of motion ( |3.8|) reduces to the massless Einstein's equation for g with an observable 
cosmological constant, 

A obs = A-6m, 2 (l-c). (3.16) 

Hence all background solutions are the same as in Einstein-Hilbert gravity with A i, s . 
However, metric fluctuations around this background are massive, with a mass term, 

\m\ {h!thl - (h^ 2 ) , ml = m 2 c . (3.17) 

The indices are raised with g. This can be read off from the minimal massive action after 
taking into account A Q b s , which is itself read off from the metric equation (not from the 
action) . 

Since c > 0, tuning it can neutralize a positive A < 6m 2 , as well as any negative A. 
Note that m 2 need not be small to recover the GR solutions and that the fluctuations 
are light as long as c is small. In the extreme case, m 2 could be of the same order as 
a quantum field theory A (say, at TeV scale) and c << 1 to cancel a large cosmological 
constant with effectively light gravitons. In the special case of A = 6m 2 , one ends up with 
a graviton mass of the same order as the observed cosmological constant, A G f, s ~ 6m|. The 
ability to cancel a A in this way generalizes the mechanism of screening of cosmological 
constant, discussed for a flat spacetime in 29, 22, p^| , to any background sourced by any 



Tuv- What is more, since the classical solutions are the same as in GR, the light bending 
discrepancy between massless and massive gravity is automatically avoided for any source, 
without invoking the Vainshtein mechanism (although this mechanism may still apply to 
the fluctuations). 

In [[34|, |35|] the authors considered massive graviton excitations over an FRW solution of 
massless general relativity. They studied the ghost related stability bounds on the Hubble 
parameter in terms of the graviton mass, generalizing earlier work on massive gravity in 



de Sitter spacetime ]3£]. It is easy to see that the ansatz (3.14) with g an FRW metric, 
provides precisely the non- linear framework in which the set up of [34, 35] is realized. Of 
course, for the reasons mentioned above, to get a consistent picture, one has to include 
appropriate dynamics for f„ v . Then the theory should provide a consistent description 



beyond the regime of classical stability breakdown [35]. 

A purpose of this discussion is to emphasize that the relevance of vDVZ discontinuity 
B P to the observational viability of massive gravity depends on the choice of f^ u . In 
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general, if this tensor is close enough to for a certain classical solution of the massive 
theory, this in turn can be very close to the corresponding solution in the massless 
theory. Thus the observational constraints imposed by the vDVZ discontinuity can be 
circumvented without invoking the Vainshtein mechanism 11 . This also underlines the need 
for further investigating the dynamics of fnu which is beyond the scope of this paper. 

4. The General non-Linear Massive Action 

In this section we show that the minimal massive action of the previous section can be 
generalized to, at most, a two-parameter family of non-linear actions. The generalizations 
are constructed in terms of a deformed determinant and correspond to adding interactions 
to the minimal theory. We elaborate on the formal structure of these terms. The generalized 
actions are ghost free at least to lowest orders in an expansion around flat space and for 
a flat auxiliary metric f^ v . Hence any ghost free theory of massive gravity must belong to 
this class. 

4.1 The generalized mass term as a deformed determinant 



The massive action ( |3,4j ) is the minimal deformation of det \/ g~ l f that satisfies the two 
criteria listed in the beginning of section 3: (i) it is free of ghosts at the lowest order in 
the fields when perturbing around flat space and (ii) it reduces to the covariant Fierz-Pauli 
mass term at lowest order in Hy . Let us now consider the most general deformation of the 
determinant with deformation parameters a n . Taking y g~ l f = K + 1 we define, 
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detv^/ = det(l + K) = £ - n ^...^a„ +1 ...a 4 c^W-A, K /* . . . K ^ . (4 .i) 

n=0 ^ 

The ordinary determinant corresponds to a n = 1. For future convenience we note that the 
deformed determinant (4.1) can be rearranged to give an equivalent expression in terms of 



\J g 1 f, rather than IK, 
l 

r!(4-r)! Cw 



r=0 



where, 



' ' 1 \ n+r 



(-1Y ■ 

/3 r = (4 - r)! u J, TT a " ' ( 4 ' 3 ) 

^ (4 — n)\(n — r)\ 

It will often be easier to manipulate the deformed determinant in this form. 

The claim is that, in terms of the deformed determinant, the full non-linear action is 
given by, 

S = -Ml j d 4 xy/=g R - 2m 2 dettvV 1 ./) + A' . (4.4) 

"Although the massive fluctuations in such backgrounds may still exhibit the vDVZ discontinuity, they 
will not affect light bending. An exception is massive gravity in AdS backgrounds which is known to not 
exhibit the vDVZ discontinuity at all |j7], 

12 For the curved space e-tensors, we use the conventions of Wald in General Relativity km, page 432. 
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This action automatically satisfies the two criteria specified above: At the lowest order in 
the metric perturbation around flat space, and for f^ v = 77^, we have IK^L „ = d^d^n. 
Then ( |4,1| ) reduces to the most general total derivative term for 7r, insuring the absence 
of ghosts at this order. Also, the generalized mass term contains the minimal mass term 
Q3.4|) for the right set of parameters and hence reduces to the covariant Fierz-Pauli action 

The five coefficients a n (or equivalently (3 n ) in fact represent only two free parameters, 
besides the graviton mass m and the cosmological constant A. The lowest order term ao 



can be absorbed into a cosmological constant as was done in (3J3). For clarity, we will 
define ao i n such a way that flat spacetime can be considered a valid background when the 
cosmological constant A is zero. This means that the mass term given by the deformed 
determinant should contain no terms linear in the metric fluctuation h. This is equivalent 
to setting ao = ot\. Then the total effective cosmological constant is given by A. 

A second coefficient can be fixed by the appropriate definition of the graviton mass 
m?. That is, we define the a n such that the term quadratic in h^ u becomes the FP mass 
term ( |2.1| ) . This is equivalent to setting a2 — ao = 1 • A third coefficient can be eliminated 
by observing that since the undeformed determinant det \J g^ 1 f does not contribute to the 
equations of motion it can be used to subtract off one term of the deformed expansion. 
We will show this below. Thus the generalized action depends on only two free parameters 
due to its anti-symmetric structure. The deformed determinant does not allow for any 
further generalization. The recursion relation constructed in 0], and the vanishing of the 
higher order terms arising from it, can be understood in terms of this structure, as will be 
explained below. 

The proof of the absence of ghosts to quartic order in non-linearities reduces to the 
calculation given in |7j and will not be repeated here. The presence of the square-root 
matrix complicates the analysis of the general case, but there are indications that the 
ghost analysis can be extended to higher orders in an ADM formalism. The general action 



in the form (4.4) also allows one to go beyond the flat auxiliary metric and consider classical 
solutions for more general f^ u . This allows for a richer class of solutions, however, now the 
absence of ghost has to be checked separately for different choices of f^ u . The structure of 
the general massive action will be discussed in more detail below. 

4.2 Formal construction 

That (^4|) is the most general massive gravity action needs a little more clarification. The 
deformed determinant contains terms that are at most fourth order in \J g -1 f . It is natural 
to ask whether the Lagrangian can contain higher order terms. These terms would have 
to vanish identically when evaluated at = in order to be consistent with the ghost- 
free construction. The epsilon tensor structure of the deformed determinant above would 
suggest that no such terms exist. However, in (?J a recursion relation is presented that might 
seem to allow one to construct higher order terms in the Lagrangian. In what follows we 
will show explicitly that these higher order terms are identically zero in 4 dimensions. 

The vanishing of these terms can be simply understood in the following way: Writing 
det(l + X) as e Trln ( 1 + x ) f or a general matrix X, and expanding the exponential in powers of 
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X gives an infinite number of terms. For a 4 x 4 matrix, only terms to order 4 are non-zero 
by the definition of determinant. These terms appear in the construction (4.1). The higher 
order terms would contribute in higher dimensions, but vanish for d = 4. These are the 
terms constructed by the recursion relation given below in ( f4.8p . 

We start by considering a general N x N matrix X. The determinant of 1 + X can be 
written in terms of the eigenvalues Aj of the matrix X, 

N N 

det(l + X) = J](l + \ l ) = Y^ e*(Ai, • • • , X N ) . (4.5) 

i=\ k=0 

The terms e&(Ai, . . . , Ajv) are elementary symmetric polynomials of the eigenvalues. That 
is, they are the sum of all distinct products of k distinct Aj, 



e (Ai, . . 


• , Aat) 


= 1, 


ei(Ai, . . 


• , Aat) 


= Ai + ... + Ajv 


e 2 (Ai, . . 


• , Aat) 


= £i<j AiAj, 


ei\r(Ai, . 


• • , Aat) 


= A1A2 ... A^ = 


e fc (Ai, . . 


• , Aat) 


= for k > N . 



(4.6) 



The eigenvalues can also be used to represent the trace of powers of the matrix X m . Here 
we will use square brackets to denote the trace, 



N 



l ]=£AT. (4.7) 



i=i 

Newton's identities give a relationship between the e^{Xi, . . . , Ajy) and the trace terms, 

k 

k 



k 

e fc (Ai, . . . , X N ) = (-!) m [^ m ] e*-m(Ai, • • • , A^) . (4.8) 



m=l 



These identities can be used to write the polynomials efc(Ai, . . . , Aat) entirely in terms of 
the traces [X m ]. For a generic 4x4 matrix one finds, 



e (X) = 1, 

ei(X) = [X], 

e 2 (X) = I([X] 2 -[X 2 ]), 

e 3 (X) = i([X] 3 -3[X][X 2 ] + 2[X 3 ]), 

e 4 (X) = ^([X] 4 - 6[X] 2 [X 2 ] + 3[X 2 ] 2 + 8[X][X 3 ] - 6[X 4 ]) , 
e fc (X) = for k > 4. 



(4.9) 



The efc(X) are exactly the terms that appear in the deformed determinant (|4.l| ) as can be 
seen from ( [4.51) . Replacing X by IK = y/ g _1 f — 1 and setting N = 4 we can write, 

4 

detv / ^ T 7=£«ne„(K). (4.10) 

n=0 
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Given this formulation, it is apparent that the recursion relation given in Q is equivalent 
to Newton's identities (4.8) with the identification C^l — > n\ e n (IK). For any N x N matrix 



this same recursion relation gives for k > N. Thus the higher order terms in y g~ l f that 
give the correct ghost free structure in the decoupling limit must be zero for n > 4. 

The existence of only two free parameters can also be seen at the level of the perturba- 
tive Lagrangian presented in ||. The first two parameters of that theory can be related to 
the a n in the deformed determinant as c 3 = 5(03 — ocq) and d§ = — jg(«4 — uq). However, 
using the above results, the terms multiplying fa can be shown to vanish. To see this, note 
that 65 (X) = for a generic 4x4 matrix. Then, from Newton's identities, the following 
combination of terms must be identically zero, 

[X] 5 - 10[X] 3 [X 2 ] + 15[X][X 2 ] 2 -20[X 2 ][X 3 ] + 20[X] 2 [X 3 ] - 30[X][X 4 ] + 24[X 5 ] = 0. (4.11) 

If we take X to be the matrix OH, then the coefficient fa in || appears in the action 
multiplying precisely this combination of terms. Analogous arguments can be used to 
eliminate free parameters that appear at higher orders in OH. Thus in four dimensions there 
are only two free physically meaningful parameters in the massive action. The higher order 
terms in an OH expansion can be determined entirely from the coefficients C3 and d*,. 

4.3 Equivalent descriptions of the action 



The full non-linear action (|4.4f) can be written in terms of the e n {K) given by (4.9) above: 



S = — M 2 f d A x^f=g~R{g) + 2M 2 m 2 I d A x^g~ ^ a n e n {K) . (4.12) 

With this formalism, it is straight-forward to convert between equivalent descriptions of 
the action. Let us first make contact with the Lagrangians presented in [p0| , B f?J. 
To eliminate the redundant parameters we rewrite, 

4 

aoeo (K) = Q det7^ T 7-^ a e n (IK). (4.13) 

n=l 

The potential term of the action becomes, 



2M 2 m 2 



(K) + a det v^/j , (4.14) 



where a n = a n — olq. As shown above, the determinant term «o det(\/ g -1 /) is non- 
dynamical and can be dropped from the action. In order for flat spacetime to be a valid 
solution, terms linear in h must vanish in the action when expanded around rjav- It follows 
that 6i\ = 0. Finally, fixing the coefficient of the mass term to be of the canonical form 
gives «2 = 1. The action ( |4.12| ) reduces to 

S = -M%J d 4 x^gR(g) + 2M*m 2 J d A x^g (e a (K) + a 3 e 3 (K) + a 4 e 4 (K)) . (4.15) 
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This action has two free parameters, 03 and 0:4. Written in this way it is easy to see that 
the action recovers the covariant Fierz-Pauli form when expanded in H . In this form it can 
be compared to the results of [^0], ^, [7) . When the parameters 0:3 and 6:4 are set to zero, 
one obtains the minimal resummed theory presented in (7j . When 03 = S4 = 1 one obtains 
the minimal deformed determinant action presented above. More generally, the parameters 
0:3 and 04 are related to the coefficients C3 and d$ of || by 03 = 6C3 and 04 = —48^. 

Because of the non-dynamical nature of the det y/ g~ x j term, the action ( |4.15| ) has an 
equivalent description in terms of the lower order e n . Moreover, we can use the definitions 
Q4.2p and (|4.3| ) to rewrite the action in terms of e n (y / g~ l f) rather than e n (K). When 
manipulating the action non-linearly it will be more convenient to work with this simpler 
description. After some rearranging, the action is given by 



S 



3 

-Ml f d 4 x^=g- R{g) + 2M 2 p m 2 f d A x^~g ^ f5 n e n { VF 1 ?) , (4-16) 



where the f3 n are given by (|4.3| ): 



f} Q = 6 — 4a 3 + a 4 , 
Pi = —3 + 3a3 — 0:4 
fh = 1 — 2«3 + «4 , 



(4.17) 



This action contains terms that are at most third order in \J g~ 1 f rather than fourth 
order as in ( 4.15| ). When expanded in h^ v this action reduces to the basic form of the 
Fierz-Pauli mass (^^) at lowest order. The minimal theory presented in the previous 
section corresponds to the choice of coefficients /3q = 3, f3\ = — 1 and 02 = 03 = (i.e., 
CX3 = «4 = 1). We stress that even though this action does not reduce to the covariant 
Fierz-Pauli form, the difference is cosmetic and the actions ( 4.15 ) and ( 4. 16| ) are physically 
equivalent (as long as is non-dynamical). 

For completeness we also present here the full non-linear equations of motion. To vary 
the action we use the relationship, 



5Tr 



■n 



Tr 



(VF T 7)'V 1 



Then it follows that, 



(4.18) 



-9 



-g< 



m=0 



Then the variation of the action gives, 

,2 3 



R 



fJ,V 



X 



n=0 



GmT. 



(4.20) 
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where Yr n \ = Y( n \(^\fg 1 /) and for X = \J 1 g 1 f we have defined (with square brackets 
again denoting the trace), 



*(0)W 



1 



x 



i[x] 



(4.21) 



V(2)(X) 



X 



■2 



X[X] + ^1([X] 2 - [X 2 ]) 



^(3)W 



X' 



•3 



X 2 [X] + ^X([X] 2 - [X 2 ]) - |l([X] 3 - 3[X][X 2 ] + 2[X 3 ]) . 



5. Discussion of results and conclusions 

Our main results are summarized in section 1. Here we discuss some of their implications, 
and related issues. 

How general is the setup? One might wonder about the generality of the non-linear 
massive gravity actions constructed via the methods of j|, ||, ||, [7j . We have argued in this 
work that in fact these methods are more powerful than they may initially appear, in more 
than one way. The actions in |, §, §, § are constructed so that they reduce to the covariant 
FP form ( [2.5D at the lowest order (quadratic in this case) in H. At first sight, they seem to 
exclude mass terms that reduce to ( |3.6|) , even though these are completely consistent with 
the basic FP mass (|2.1| ). However, we have shown that the excluded actions differ from 
the rest only by non-dynamical terms and, in this sense, are not left out. Also, within our 
formulation the origin of the 2-parameter family of actions, and the impossibility of adding 
more parameters, becomes manifest. 

Moreover, mass terms considered in ^ ^] and here, depend only on scalar functions 
of the form F{g~ l /). One may wonder about mass terms that involve scalar functions 
F'(f~ 1 g). As pointed out in section 2, the quintic order ghost analysis of || is general 
enough to allow for both forms. However, it so happens that at the non-linear level, 
F(g _1 /) is adequate to express all admissible (by the no-ghost condition) non-linear mass 
terms in a simple way. 

Of course, the theory presented here can be extended by introducing dynamics for /, 
considering additional auxiliary metrics or other such modifications. However, any theory 
of massive gravity that reduces to the basic FP action at lowest order in h, is ghost-free to 
lowest order in interaction, and depends only on the metric g and the auxiliary metric /, 
must belong to the two-parameter family of actions presented here. 

Non-flat auxiliary metric. The ghost analysis of massive gravity actions is valid 
only for a flat auxiliary metric. The associated constraints introduce deviations from 
general relativity. For such theory to comply with observations, at the least one needs the 
Vainshtein mechanism On the other hand, if one allows the auxiliary metric 

to be non-flat, or better, related to g^ as in section 3.3, then the classical solutions could 
be very close to general relativity while still retaining some attractive features of massive 
gravity, like the screening of the cosmological constant. The proper arena of realizing such 
solutions is theories with dynamical f^ u , though not necessarily restricted to bi- gravity 
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theories. The ghost issue in such models has to be re-examined. In fact, the recent analysis 
of massive gravity in FRW spacetimes [34, ^] is relevant to this construction. 



The ghost issue. It is clear that the appearance of the square-root matrix \/ g~ 1 f 
is an inevitable feature of the Goldstone boson formulation of massive gravity Q. Such a 
matrix cannot be easily manipulated analytically. Thus the very setup that simplifies the 
ghost analysis in the decoupling limit, introduces a complication at the non-linear level. 
A straightforward ADM analysis of the ghost is therefore not an option, although such an 
analysis can be carried out for some special cases, for example, for JVj = 0. For more on 
this, see @. (also see the note added below) 

Another issue is that even for models with flat f^ v , where the ghost analysis applies, 
the 2-parameter family of actions may need to be further constrained to a single parameter 
if one is to avoid unstable non-asymptotically flat solutions, as can be seen from [29]. Thus, 
while the candidate ghost-free actions considered in || [7| and here are a step forward, the 
final question of non-linear massive gravity is far from settled. 

Note added: After this work was submitted to the arXiv, [[5]] appeared with some related 



results. More recently the absence of ghost at the non- linear level was shown in 41]. 
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